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The so-called V15 molecule, of formula Ke[Vi5 Ase042(H20)]8H20, is formed of 15 spins 
i antiferromagnetically coupled. The resultant spin is equal to i. Contrary to what is 
expected at first sight, this half-integer spin is gapped. We show that this is a consequence 
of the multi-spin character of the molecule. Time reversal symmetry (Kramers theorem) is 
preserved due to the existence of a groundstate degeneracy larger than 2. Magnetization ex- 
periments have been performed at different temperatures, sweeping-field rates and couplings 
with the cryostat. Their interpretation in terms of the Landau-Zener-Stiickelberg model 
gives evidence for adiabatic LZS spin reversal with and without dissipation. 

§1. Introduction 

Contrary to high spin molecules such as Mni2 or Fe§ (S = 10) 1 \ the V15 
molecule has a low resultant spin (S = ^). This leads to important differences 
with respect to high spin molecules. In particular, energy barrier between spin-up 
and spin-down directions, as well as dipolar couplings between molecules, become 
negligibly small. Nevertheless, in low spin molecules such as V15, the Hilbert space 
is as huge as in large spin molecules. This is simply because they are both multi- 
spins. In V15 the Hilbert space dimension is equal to 2 15 ; in Mni2 and Fes it is 10 8 
and 6 8 respectively. In principle, this molecule should be discussed in terms of the 
entanglement of the 15 spins ^. However, considerable simplifications occur at low 
temperature. The groundstate spin is equal to \ and the first excited one is equal to 
|. These levels being reasonably well separated in energy, low temperature experi- 
ments deal with the groundstate only. We showed experimentally that in zero field, 
the symmetric and antisymmetric superposition of the groundstates | |, S z = i) and 
S z = — 5) are split by an energy A. The origin of this splitting with a half-integer 
spin is attributed to its multi-spin origin. In the absence of decoherence, and in zero 
field, the ensemble of molecules should coherently oscillate between the symmetric 
and antisymmetric states | ^, S z = \) ± ||, S z = —\) with the frequency A/h. In the 
presence of a sweeping field vh = dH/dt, the initial state \^,S Z = ^) either rotates 
to the state |^, S z = — ^) or is not modified, depending on the sweeping field rate. 
The switching probability depends on the comparison of the time r sw = A/vng/J-B 
during which the field sweeps the mixing region (~ A/gfi^), with the oscillation 
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time t os = h/A. When the field changes slowly (vh <C A 2 /figns) the system stays 
in the same state (groundstate) and the spin rotates on an adiabatic way from 3) 
to 1 1,— |). When the field changes rapidly (vh 3> A 2 /H/j,b) the system has not the 
time to experience quantum mixing and the spin does not rotate. The switching 
field probability is given by P = 1 - exp(-r) with r = irA 2 /2hv H gHB- This is the 
main result of the Landau-Zener-Stiickelberg (LZS) without dissipation. This model 
was introduced for inter-band tunneling in semiconductors 2 \ It has been recently 
applied to the specific case of mesoscopic magnetism by Miyashita 3 ). The effects of 
dissipation have been considered by several authors 4 ). 

Our time resolved magnetization measurements performed at different conditions 
of field, temperature and coupling with the cryostat, are interpreted on the basis of 
a LZS model. Following Abragam and Bleaney, we introduced the role of phonons 
and spin bath through simple spin-phonon and spin-spin transitions 5 ) . 

After this introduction, the paper is ordered as follows: 2) Structural and ther- 
modynamic magnetic properties, 3) Time resolved magnetization measurements, 4) 
Quantitative approach and evidence for gapped of half-integer spin, 5) Origin of the 
gap in the V15 molecule, 6) Conclusion. 

§2. Structural and thermodynamic magnetic properties 

The so-called V15 molecule is a complex of formula K 6 [V I 1 YAs 6 4 2(H20)]8H 2 0. 
This complex forms a lattice with trigonal symmetry (a = 14, 02.9 A, a = 79.26°, 
V = 2632 A 3 ), containing two V15 complexes per cell. The third order symmetry 
axis of the unit-cell is also the symmetry axis of the V15 clusters (space group R3c). 
All the fifteen V IV ions of spin S = \ are placed in a quasi-spherical layered structure 
formed of a triangle sandwiched between two non-planar hexagons (Fig.l). 

There are five antiferromagnetic exchange constants. Each hexagon contains 
three pairs of strongly coupled spins (J= 800 K) and each spin at a corner of the 
inner triangle is coupled to two of those pairs (one belonging to the upper hexagon 
and one belonging to the lower hexagon). We have two different determinations of 
the exchange couplings shown in Fig.l: J\ = J' = — 30 K, J 2 = J" = — 180 K, 
J = -800 K 6 ) or Ji = J' = -150 K, J 2 = J" = -300 K, J = -800 K 7 \ These 
values should not be taken as definitive. There are too many free parameters in these 
evaluations. However, the orders of magnitude are correct. The V15 molecule can 
be seen as formed by three groups of five V IV ions. Each group with a resultant spin 
^ is located on a corner of the inner triangle. These three spins 1 ; interact with each 
other through two main paths, one passing by the upper hexagon and one passing 
by the lower one. They are related to each other by the threefold symmetry axis of 
the molecule. This is a typical example of a frustrated molecule. 

The magnetization measurements given below were obtained on single crystals 
of the V15 complex. A small dilution refrigerator allowing measurements above 0.1 K 
was inserted in an extraction magnetometer providing fields up to 16 T, with low 
sweeping rates of about 1 minute/point. Below 0.9 K we clearly observed three 
jumps, one in zero field, and two at H\ = ±2.8 T (Fig.2). They correspond to the 
— \ 44> \ and ±.\ to ±| spin transitions, with respective saturation at 0.50 ±0.02 \xb 
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and 2.95 ±0.02 per V15 molecule. Above 1 K, the observed jumps vanish because 
the S = I level becomes thermally occupied. The magnetization curves are reversible 
and do not show anisotropy when the field is applied along different directions, within 
the accuracy of the experiments. Furthermore, experiments made on a single crystal 
or on several non-oriented single crystals were the same. So, any energy barrier 
preventing spin reversals must be quite small (less than 50 mK). The magnetization 
curves measured at equilibrium (Fig. 2) are fitted using the Heisenberg Hamiltonian 
for three frustrated spins S = \ : 

H = — ^2 J a (S a lSa2 + S a 2S a 3 + S a 3S a i) — g[lBH(Si + S2 + S3) (1) 

a=x,y,z 

where Si = S2 = S3 = J a < and g = 2. 

The three spins \ represent the resultant spin of the three groups of five spins 
(Fig.l). The energy scheme for the lowest levels is obtained by diagonalization of the 
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Fig. 1. Quasi-spherical layered structure formed 
of a triangle of V IV sandwiched between 
two non-planar hexagons. Main exchange 
paths have been indicated. Scheme of the 
5 spins model used for our calculations. 
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Fig. 2. Magnetization curves measured on a 
single crystal of V 15 . The agreement with 
the calculated curves obtained with our 5- 
spins model is excellent for Jo = — 2.445K 
(continuous curves). However, the spin 
transition measured at 2.8 T and at low 
temperature is broader than the calculated 
one. Anisotropic exchange interactions do 
not influence the shape of the transition 
(dashed curve for J x = J y — — 2.75K, J z = 
— 2.2K). This could be done by antisym- 
metric D-M interactions. 



3-spins matrix of the Hamiltonian 1. This matrix contains 64 terms, and the chosen 

basis is: {| + ++), | ++-), | + -+), | -++), | + ), | -+-), | - -+), | )}. In the 

isotropic case and in zero field, the first excited level with spin | has, as expected, 
a fourfold degeneracy. However, it turns out that the ground state, with spin ^, 
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has also a fourfold degeneracy. This is a consequence of its multi-spin character: 
the groundstate of V15 is made of the juxtaposition of two independent spins 5, 
each with a degeneracy of 2. We will see below that this "degeneracy doubling" 
effect has a general consequence on the quantum reversal of a spin ^ and more 
generally of a half-integer spin, if interactions are antisymmetric. It is easy to see 
analytically that the separation between the ground state and the first excited state 
is —3 Jo/2, and that the value of the field at which the configuration S = Sz = § 
becomes favorable against the S = Sz = \ is U\ = — 3 Jq/ {2g^). The measured 
field Hi = 2.8 T gives Jo ~ —2.51 K and then a separation of « 3.76 K. These 
results are confirmed by high temperature susceptibility experiments (see below). 
The value of « 3.76 K is extremely small, compared to the exchange interactions 
Ji,J2, J- This is a consequence of the spin frustration of the molecule. The three 
spins ^ of the sandwiched triangle are coupled via the two hexagons, as indicated in 
Fig.l. According to this scheme, the exchange fields at the site of each triangle spin 
should cancel each other. The observed small value of Jo results from the exchange 
couplings Ji and J2 not cancelling each other exactly. 

This is why the separation between the lowest levels S = \ and S = § is al- 
most two orders of magnitude smaller than these interactions. The magnetization 
is calculated from the projection of the magnetic moments upon the applied field 
Szi = {i\Szi + Sz2 + Sz3 \i) , averaged over the different eigenvectors \i) of energies 
Ei. The population of different energy levels is taken at equilibrium (Boltzmann 
distribution). In the isotropic case, an analytical form of the magnetization is ob- 
tained 8 ^. The fit of the magnetization curves measured at different temperatures 
gives Jo ~ — 2.445K, which is very close to the value of 2.51 K independently ob- 
tained above. The agreement between calculated and measured curves is excellent at 
each temperature, except for the slope of the transition at 2.8 T which is broader in 
experiments. This broadening of about 0.7 T, does not depend on samples and can- 
not be ascribed to dipolar or hyperfine field distributions (about 1 mT and 40 mT 
respectively), nor to anisotropy exchange effects included in some of our calcula- 
tions. However, we found that antisymmetric Dzyaloshinsky-Moriya interactions 
Ddm {Si x Sj) (Ref. 9 and refs. therein), which are allowed by symmetry and 
couple the states S = \ and S = |, might explain this broadening (see section 5). 
In fact the important effect of Dzyaloshinsky-Moriya interactions for the induction 
of tunnel splittings, otherwise forbidden, was pointed out in 10 ) for the case of the 
Mni2 molecule. Detailled calculations have been done in n ) . 

The susceptibility is calculated by taking the first derivative of the analytical 
expression of the magnetization 8 ) . The zero-field susceptibility \i has a very simple 
expression: 

{gUB? 5 + e-^ kT 
Xl Ak B T 1 + e -3Jo/2fcT W 

The limits T <C 3Jo/2/cs and T 3> 3Jo/2/cb correspond respectively to a spin ^ with 
the effective moment \/3<7/zb/2, and 3 spins \ with the effective moment 3(7^/2. The 
fit of the susceptibility curves measured between 0.1 K and 100 K is given Fig.3 and 
Fig.4. The measured reciprocal susceptibility vs. temperature and its comparison 
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with Eq.(2) for Jo ~ —2.445 K (determined from M{H)) is excellent. The straight 
line represents the Curie- Weiss law with the Curie constant C = 0.686 /xbK/T and 
the paramagnetic temperature 6 = —12 mK. This Curie constant corresponds to 
the effective paramagnetic moment 1.75 /xb, a value very close to gn^(S(S + l))i = 
<7/xb\/3/2 for S = 5. The paramagnetic Curie temperature is of the order of dipolar 
couplings between S = \ molecule spins. Above 0.5 K, the susceptibility deviates 
from the experimental data and from the calculation from Eq.2. This is because the 
spin ~ ground-state starts to be deteriorated. Fig.4 shows the measured effective 
moment fj, e g = ^3kT\i as a function of temperature. It increases from the value 
of 1.75 /Ub corresponding to the spin | ground state, and tends asymptotically to 
the value of 3 /iB which is preserved up to 100 K (three independent spins ^). The 
fit of these susceptibility measured between 0.1 and 125 K to Eq.2 is excellent with 
Jo = -2.445 K. 
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Fig. 3. Measured and calculated reciprocal 
susceptibility vs. temperature. The dashed 
curve comes from the model with Jo = 
— 2.445 K and the continuous line is a 
Curie- Weiss fit with a paramagnetic Curie 
temperature 9 = 12 mK and a Curie con- 
stant corresponding to an effective moment 
of 1.75 fi B (S = I). 
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Fig. 4. Curve calculated in the 5-spins model 
is valid up to 125 K, showing that excited 
levels corresponding to triangle-hexagons 
couplings are well separated from the trian- 
gle levels. Two plateaus can be seen: one 
at 1.75 fiB corresponding to one spin | (be- 
low 0.5 K) and one at 3 /j,b corresponding 
to three spins \ (above 25 K). 



§3. Study of the two-level system realized in the V15 molecule 

In this section we concentrate on the zero-field \^,S Z = \) 44> \^,S Z = —^} 
transition. We studied the temperature and time dependence of a small single- 
crystalline grain of V15 (L ~ 30 — 40 /*m), with different couplings to the cryostat. A 
first comparison with the LZS model shows that with zero-field splitting of the order 
of 50-200 mK (see section 4) the transition is purely adiabatic if vh "C 10 10 T/s. 
In our experiments where vh < 0.14 T/s (see below) it is a fortiori the case. Our 
experiments being done at finite temperatures, we found a much richer behavior than 
simple adiabatic LZS, with two limiting cases in which the magnetization curves 
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M(H) are reversible, and an intermediate one where it is irreversible. This behavior 
results from effects of the environment on the adiabatic LZS model. 

3.1. Fast sweeping rate / weak coupling to the cryostat, and adiabatic LZS. 

The first limit of fast field variations and weak coupling to the cryostat, was 
reached by inserting a plastic sheet, much thicker than the sample, between the 
sample and the sample-holder and restricting the temperature to lowest values. 
The four M(H) curves of Fig. 5 show either a weak hysteresis or no hysteresis at 
all, depending on the competition of temperature (giving hysteresis through spin- 
phonons transitions) and sweeping field rate vh = dH/dt (preventing these tran- 
sitions). This effect can be accounted for by the ratio avn/T 2 (see the model be- 
low). For avn/T 2 < 14 T/sK 2 weak hysteresis is observed on the M(H) curve. 
It is attributed to finite but small probability for spin-phonon transitions. For 
avfj/T 2 > 14 — 40 T/sK 2 the M(H) curve becomes reversible. Spin-phonons 
transitions have not the time to take place and the probability for their occur- 
rence is negligibly small; only the groundstate is occupied. The \\,S Z = \) to 
\\,S Z = —\) LZS transition is adiabatic. The adiabatic M(H) curve is nicely 

fitted with the expression M = {\){g^ B ) 2 H/[{A ) 2 + (gv B H) 2 ]i obtained from 

E = ^dAn/dH where Ah = [(g^H) 2 + (^o)l^ ^ s tne two-level field-dependent 
splitting energy. The saturation (^)gfi-B is very close to 1)Ub as expected, and the 
initial slope M/H = {g^) 2 /2Aq is proportional to 1/Aq, giving immediately an or- 
der of magnitude for the zero-field splitting Aq « 0.1 K. The fit of the whole M(H) 
curve, as well as the 1/H 2 approach to saturation gives Aq ~ 80 mK. 

In short, a single crystal of V15 molecules with reduced couplings to the cryo- 
stat, is a realization of the adiabatic LZS model applied to magnetism 3 ). The spin 
temperature can then reach very low values T s = Tj/(1 + (gfiBH/Ao) 2 ) « Tj/50, 
where Tj is the initial temperature. The same procedure applied to the 5 to | spin 
transition which occurs in a rather high field (2.8 T) should allow to obtain high 
nuclear spin polarisation. 

3.2. Low sweeping rate / strong coupling to the cryostat, and dissipative LZS. 

Here the coupling between the sample and the cryostat is made better in or- 
der to favor spin-phonon transitions. We considered two cases. In the first one, 
the coupling was obtained by simple contact using some grease between sample and 
sample holder. In the second case the coupling was better, the contact being made 
of a mixture of silver powder with grease. In both cases non-equilibrium hysteresis 
loops are observed at fast-sweeping rate, down to the lowest measured temperature 
of 50 mK. The M(H) curves obtained with the first type of coupling (grease con- 
tact), can be found in Ref. 8, 9, 13. Those obtained with the best coupling with 
the cryostat (Ag+grease) are similar. They are shown Fig. 6. When the field in- 
creases, coming from negative values, the magnetization passes through the origin 
of the coordinates, reaches a plateau and then approaches saturation. This leads to 
a winged hysteresis loop characterized by a reversible but temperature and sweep- 
ing field-dependent zero-field susceptibility. The irreversible wings depend also on 
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5. Hysteresis loops measured in the case 
of a sample rather well isolated from the 
cryostat (a = 130 sK 2 ). The curve at 
60 mK and 0.14 T/s shows an adiabatic 
LZS spin transition. It is very well fitted 
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Fig. 6. Hysteresis loops measured in the case 
of a sample rather well coupled to the cryo- 
stat (a = 0.09 sK 2 ). 



with the expression 
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temperature and sweeping field. The sweeping field dependence of the zero-field 
susceptibility is typical of the phenomenon studied here. 

In Fig. 6 the hysteresis of the M(H) curves disappears when the field is slow 
enough. This is a natural consequence of thermal equilibrium. Spin-phonon transi- 
tions are sufficiently probable so that the two-level populations equilibrate with the 
Boltzmann distribution. As an example, equilibrium is almost reached at 60 mK and 
vh < 0.14 mT/s (approximated by the median of the two branches of the M(H) 
hysteresis loop, Fig. 6). This reversibility is very different from the one discussed in 
3.1 where field variations were fast enough to suppress spin-phonons transitions com- 
pletely, giving a highly non-equilibrium groundstate at finite temperature (adiabatic 
LZS). 

In order to give a first qualitative analysis of these "butterfly hysteresis loops" 
9 ) we will see how the level occupation numbers vary when sweeping the external 
field (Fig. 7). The spin temperature Ts is such that nxjn^ = exp(Ajj /&b2s) where 
"1,2 ("i,2equ) is the out-of-equilibrium (equilibrium) level occupation numbers. The 
magnetization curves at 60 mK (Fig. 6), show a spin temperature significantly lower 
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than the bath temperature (ni > ni equ , T$ < T) between 0.3 T and 0.15 T (the 
field at which the magnetization curve intersects the equilibrium one). After this 
intercept, Ts is larger than the bath temperature (ni < rei equ , T$ > T) and at 
sufficiently high fields (about 0.5 T) it reaches its equilibrium value (ni = ni eq u, 
Ts = T). Note that the magnetization curves measured between 0.7 and 0.02 T at 
fast sweeping rates are nearly the same, suggesting weak exchange with the bath, i. e. 
nearly adiabatic demagnetization. This is because below 0.5 K, the heat capacity of 
phonons C p h is very much smaller than that of the spins C5, so that the phonons 
temperature T p h very rapidly adjusts to that of the spins Tg. The spins and the 
phonons can be seen as a single system. The coupling of this system with the 
cryostat is weak because energy transfers through phonons must obey the drastic 
condition fiuj = Ah and the number of phonons available at this energy fiuj is 
extremely small: ny = / a(u)du/ exp(huj/kT). The integral is taken over the level 
linewidth Auj (due to fast hyperfine fluctuations), a(u)duj = 3yw 2 du;/27rVp h is the 
number of phonon modes between u and u + du> per molecule of volume V, and v p h 
is the phonon velocity. Taking the typical values for Vis-like systems v = 3000m/s, 
and Auj = 10 2 MHz, we find at T=0.1 K, n T « 10~ 6 - 10~ 8 phonons/molecule. 
The number of such lattice modes being much smaller than the number of spins, 
energy transfer between phonons from the cryostat and spins must be very difficult, a 
phenomenon known as the phonon bottleneck 5 )' 12 \ Despite a good thermal contact 
between the sample and the cryostat, the energy flow from the latter is not sufficient 
to compensate the lack of phonons at energies h(ui ± Aoj). Most molecules are out 
of equilibrium, at a temperature Ts different from the cryostat temperature T. 



E 




Fig. 7. Schematic representation of the V 15 spin | energy scheme in a static magnetic field (the 
level broadening is due to nuclear spin fluctuations). Probability transitions P in the LZS 
model and spin-phonon transitions are indicated by different types of arrows. The existence of 
a zero-field splitting Aq with a spin i is related to the fact that in V 15 , this scheme is doubly 
degenerated with the superposition of two independent and identical schemes. 
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§4. Quantitative approach and numerical calculations 

In a spin \ two-level system, the populations of spins and phonons obey a set 
of two differential balance equations (e.g., see Ref. 5 )). As discussed above, at low 
enough temperature, the condition b = C s /C p h 3> 1 implies T p h — ► T s / T '. The 
phonon population is dragged by that of the spins, and the problem reduces to a 
single equation. Writing x = (n\ — «2)/(ni eq — n eq ), the spin population is given 
by Tndx/dt = 1/x — 1. The solution of this simple equation is —t/rn = x — xq + 
ln[(x — l)/(xo — 1)] where th = &T p h = (a/A 2 H ) tanh 2 (Z\#/2/cBT) is the relaxation 
time of the spin-phonon system to the cryostat temperature and xq the initial spin 
population. The constant a is formally proportional to the size L of the sample, 
i.e. to the distance the phonons have to travel to reach the cryostat temperature 9 ^. 
In experiments, it is common to change L artificially by modifying the thermal 
coupling at the sample/sample-holder interface (see above). The above expressions 
of x and tjj allow to calculate the magnetic relaxation and the magnetization curves 
M(H,T,t) with the use of only two parameters Aq and a. Note that Tndx/dt could 
be written THVH(dx/H), suggesting that the parameter thvh oc avn/T 2 is relevant 
for M(H) curves. 

The M(H) curves have been calculated to give the best agreement with the 
measured ones. The obtained parameters are Aq ~ 50 mK and a ~ 0.15 sK 2 (grease 
contact). The comparison between the magnetic relaxation curves measured in the 
same conditions and calculated (not shown) gives Aq ps 150 mK. In the case of the 
best coupling between sample and cryostat , the agreement between measured and 
calculated M(H) curves gives Aq 80 mK and a ps 0.09 sK 2 . This smaller value 
of a is consistent with a better contact to the cryostat. These determinations of Aq 
must be compared to the value Aq ps 80 mK obtained from the fit of M(H) without 
dissipation (section 3.1). Relaxation time measurements being less accurate than 
M(H) fits, we can consider that the zero-field splitting of V15 is Aq ps 80 ± 20 mK. 

Although these quantitative deter- 
minations of Aq are not direct (neutron 
scattering experiments are in progress, 
see also Ref. 14), the obtained finite 
value is qualitatively proved by the 
sweeping field dependence of the initial 
susceptibility and in particular by the 
clear observation of an adiabatic LZS 
regime, in fast fields. Quantitatively, 
the value of about 80 mK seems rather 
safe for the following reasons. We ob- 
tained three couples (Aq,ol) indepen- 
dently from the comparison of measured and calculated M(H) curves in the adia- 
batic LZS regime without dissipation (bad coupling with the cryostat, a ~ 130 sK 2 ) 
and with dissipation (good and very good couplings, a ps 0.15 and 0.09 sK 2 ). The 
Aq values obtained are close to each other and the values are consistent with the 



Table I. Relaxation times at different fields 
and temperatures, obtained from the fit 
of the curves shown Fig. 8 (tj/cx P ) and cal- 
culated from the time dependent equation 
given in the text with a — 130sK 2 and 
A = 80mK. 



T(K) 
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TH cxp(s) 


TH calc(s) 


0.05 
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8716 
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14 


551 
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70 


3883 


3675 


0.15 


70 


970 


997 



10 B. Barbara, I. Chiorescu, W. Wernsdorfer, H. Bogge, and A. Muller 



whole physical picture, giving credit to our approach. Furthermore, our model is 
self-consistent. In the case of a ~ 0.15 sK 2 we have found a relaxation time of 
about 1 s, from the above expression of tjj, which is close to the measured one 13 ). 
Moreover, we obtained (i) a sound velocity v p h 3000m/s, which corresponds to 
the values given for insulating systems % and (ii) the calculated specific heat ratio 
b = THV p h/L 10 8 is comparableto the expected value for an ensemble of spins 
\ in a low temperature bottleneck regime 5 )' 12 - 1 . Finally, in the other limit where 
a ~ 130 sK 2 and Aq ~ 80 mK, our magnetic relaxation measurements are well fitted 
to the model given above and the obtained relaxation times are very close to the 
calculated ones for the whole range of fields, excepted near zero field where the mea- 
sured relaxation (Fig. 8) is about six times faster (Table I). This shows that when 
the relaxation takes place within the zero-field splitting Aq, the phonons bath is 
no longer sufficient to explain the dissipation. An excess relaxation probably comes 
from the Dzyaloshinsky-Moriya interactions and the spin bath, in particular from 
the fast fluctuations of the 51 V nuclear spins. A more detailed study of the inter- 
play between nuclear and electronic spins, in the presence of the spin mixing of the 
Dzyaloshinsky-Moriya interactions should be done in the case of this molecule, in 
the light of our recent works on Hq + ions 15 ) . For a discussion see the next section. 
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Fig. 8. Relaxation of an isolated sample, in a field H <C Ao for which spin up and spin down 
states are mixed (left), and in a field H 3> Ao where the level structure is dominated by the 
Zeeman effect. The fits (continuous curves) of data (dots) are good except when the field and 
the temperature are small (0.15 K, 0.014 T). 



§5. Origin of the gap in the half-integer ~ spin of the Vis molecule 

It is well known that time reversal symmetry prevents half-integer spins to be 
gapped in zero field (Kramers theorem 16 )). This is also the reason for the absence 
of Haldane gap in magnetic chains with half-integer spins 17 ). Our first evidence 
for a gap in spins i V15 molecules is at first sight in contradiction with this first 
principle. The difference between the spin I of V15 and an atomic spin ^ is obviously 



The V15 Molecule, a Multi-spin Two-level System: 



11 



the multi-spin character of V15 . The groundstate spin % is the resultant of exchange 
interactions as shown section 2. As we suggested in Ref. 9 and references therein, 
antisymmetric Dzyaloshinsky-Moriya (D-M) interactions can provide a finite gap in 
a multi-spin system with resultant half-integer spin. On the point view of symmetry, 
the molecule V15 is a good candidate to show this effect, due to the lack of inversion 
center. This is particularly true for spin pairs situated near the borders of the 
molecules. The Hamiltonian 1 should therefore contain the additional term: 

Hdm = —^^DijSi x Sj, (3) 

ij 

the summation being done over all spin pairs. The norm of the vector D can be eval- 
uated from the measured anisotropy of the g-factor. Very recent EPR experiments by 
Ajiro et al. 18 ) made accurate measurements of the angular variations of this factor, 
in a plane containing the c-axis. The results, much more complete than the previous 
ones 6),7) , give at 2.4 K, g\\ c = 1.98 and g ±c = 1.968. Then using D « JAg/g, with 
Ag/g = (g\\ c — 9±c)/g = 0.6%, we can see that the value DS 2 D/A m 80mK 
derived from measured Aq, gives J ~ 60K, which is of the right order of magnitude 
for average exchange interactions. Therefore D-M interactions could really be at the 
origin of the observed zero-field splitting. 

The best theoretical proof for that was given recently by Miyashita and Na- 
gaosa 19 ) who calculated the magnetization curves and splittings resulting from level 
mixings in the 3-spins Heisenberg model (Hamiltonian 1 in section 2) with various 
types of perturbations of the form 

Hi = -^ctijSizSjx. (4) 

They showed that the groundstate is not gapped, as expected from the Kramers 
theorem 16 -*, only if the aij show the reflection symmetry, i.e. otij = ctji. In other 
cases, when otij / ctji the perturbation in Eq.4 creates a gap between the ±1/2 spin 
states, as observed experimentally. However, as discussed above, the spin | being a 
multi-spin system, its groundstate degeneracy is larger than 2. In the present case 
of 3 spins \ it is equal to 4. The four states with S = \ form two sets of avoided 
level crossing at zero field, each with two-fold degeneracy. This allows Kramers 
time reversal symmetry to be preserved even in the presence of a zero-field splitting. 
Coming back to Eq.3, one can write H y = —Dij y (Si x Sj z — Si z Sj x ), which identifies 
to Hi in Eq.4 if ccjj = — ctji = ^ijy, which is not surprising as D-M interactions are 
antisymmetric. 

Beyond the case of the V15 molecule, D-M interactions, which are very gen- 
eral with molecules (and in particular if they are of finite size), should lead to a 
gapped groundstate in non-integer multi-spin molecules. It is possible that other 
mechanisms, such as the multi-exchange couplings in spin liquids, could lead to sim- 
ilar effects. The molecule V15 can be considered as an example of nanometer scale 
spin- liquid. 

Before ending this section, let us consider some characteristics 20 ) of the nu- 
clear spins of 51 V. The natural abundance of this element with nuclear spin |, is 
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99.76%. The resonance frequency and magnetogyric ratio of 2.63 GHz and 7.045 
rad/10 7 Ts give an hyperfine field of about 40 T. The inverse hyperflne field act- 
ing on electronic spins, roughly 2000 times smaller, is in the 20 mK scale, i.e. of 
the order of Aq ~ 80 mK. As indicated before, this field will give a broadening 
of electronic levels. This broadening corresponds to the incoherent bunching of 8 
electro-nuclear levels. Another origin for level broadening comes from the spin mix- 
ing by Dzyaloshinsky-Moriya interactions between the different pairs of spins in the 
molecule. The observation of faster relaxation near zero field, suggests that fluc- 
tuations on the two levels intercept each other near zero field (at the energy scale 
of Aq). This is the case for nuclear spin fluctuations, the broadening of which is 
of the order Aq/2. These fluctuations should therefore be relevant to explain the 
excess of relaxation in zero field 8 )> 9 ). Level broadening, resulting from incoherent 
nuclear spins or Dzyaloshinsky-Moriya fluctuations, can be considered as a magnetic 
external noise on the LZS model 21 ). Such a noise cannot explain the origin of Aq, in 
particular because it cannot, by itself, simulate the LZS model, nor the associated 
dynamics in the presence of spin-phonon transitions. As shown above the origin of 
Aq is clearly associated with the static antisymmetric effect of Dzyaloshinsky-Moriya 
interactions. 

§6. Conclusion. 

The V15 molecule, a multi-spin system with spin ^, shows clearly adiabatic LZS 
transitions with or without dissipation, depending on the sweeping field rate and on 
the coupling with the cryostat. This allows to study environmental effects on quan- 
tum spin reversal in a two level system. In particular it is shown that the observed 
magnetic relaxation is associated with spin-phonon transitions in the Zeeman regime 
(phonon bath). Near zero field, level broadening due to spin fluctuations (essentially 
nuclear spin fluctuations) seem to overlap, giving an excess of relaxation, associated 
with the spin bath. An important consequence of this study is that, contrary to 
what is expected at first sight, this half-integer spin is gapped in zero-field. This is a 
consequence of the multi-spin character of the molecule plus Dzyaloshinsky-Moriya 
interactions. Time reversal symmetry is preserved due to a fourfold degeneracy of 
the groundstate. 
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